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■ Abstract 

We use higher derivative classical gravity to study the nonlinear coupling 
between the cosmological expansion of the universe and metric oscillations of 
Planck frequency and very small amplitude. We derive field equations at high 
orders in the derivative expansion and find that the nature of the new dynamics 
{Sj | is extremely restricted. For the equation of state parameter w > the rela- 

tive importance of the oscillations grows logarithmically. Their effect on the 

■ cosmological expansion resembles that of dark energy. 

o 

m 

^ '. 1 Introduction 

p 

A general gravitational action can be arranged as a derivative expansion, where higher 
derivative terms contain more factors of the curvature tensor and its covariant deriva- 
tives. As such this could be viewed as the expansion of some underlying theory. But 
major problems arise if some truncation of the derivative expansion is treated as a 
fundamental theory. For example the inverse propagator of the high derivative theory 
is at best of the form 

k 2 Y[(k 2 - m 2 ) with all m 2 > 0. (1) 

i 

Close to every other zero the k 2 will have the wrong sign. These ghost poles lead 
to loss of unitarity, the possibility of unstable solutions, etc. At a deeper level the 
rrii probably have little to do with the spectrum of the underlying theory. Also 
a truncated derivative expansion will have solutions bearing no resemblance to the 
solutions of the true theory. Higher derivative gravity only makes sense as an effective 
field theory where the higher derivative terms are properly treated as perturbations 
|T||. (For a contrary view see ]3[].) 
*bob.holdom@utoronto.ca 



i 



Nevertheless we will show in this paper that higher derivative gravity in the case 
m% > can serve as an instructive and suggestive framework for the investigation of 
the following question. Can there be a dynamical nonlinear coupling between physics 
at very large and very small time scales (UV-IR coupling) in a theory of gravity? 

We will consider a class of solutions to higher derivative gravity theories of the 
FRW form, where the standard FRW time dependence for the scale factor can be 
recovered as a special case. This special case corresponds to particular initial con- 
ditions, and departures from this will yield oscillations with Planck scale frequencies 
in the scale factor of the FRW metric. It has been noted before [0] that rapid met- 
ric oscillations can emerge from higher derivative theories. We will find that the 
oscillations, even when averaged over, produce a departure from the FRW solution, 
where the fractional amount of this departure grows with time. We refer to this as 
a weak dynamical instability since the growth of these oscillations turns out to be 
logarithmic. 

We claim that this is a generic property of higher derivative gravity, true at any 
order in the derivative expansion. We have reached this conclusion through a study 
of the field equations at high order in the derivative expansion. Since Planck frequen- 
cies are involved, derivatives of the oscillating component are not suppressed. The 
amplitude of the oscillation on the other hand is extremely small, at least at late 
times, if the impact of the oscillations on the equations is smaller than or of order 
the usual effects of cosmological expansion. We therefore consider a small amplitude 
expansion, where terms quadratic in the scale factor and its derivatives produce the 
leading nonlinear effects. It is these nonlinear effects that couple the oscillations and 
cosmological expansion together and give rise to the dynamical instability. 

As we have said the new solutions can be arbitrarily close to the standard FRW 
solution, in which case there will be an arbitrarily long time where the effects of 
all higher derivative terms are much smaller than the Einstein curvature term. But 
eventually the effects of higher powers of curvature can become as important as the 
Einstein term, and the effect on the cosmological expansion can be significant. 

One would naively think that the large and growing number of parameters in 
higher derivative actions would produce results that are very ambiguous. On the 
contrary we find that in the small amplitude expansion that the derivative expansion 
is extremely restrictive. Each order in the derivative expansion, at six derivatives and 
beyond, results in only two parameters, one of which multiplies both the leading and 
next-to-leading terms in the small amplitude expansion. And then quite surprisingly 
we find that the second parameter drops out of the quantities of interest. Thus a 
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very particular pattern emerges in the derivative expansion that allows us to relate 
the amplitudes and energies of the oscillating modes to the matter sources. 
We start with the gravitational action 

S = M%iJ d 4 xy/=g~(-2A + R + aR 2 + bR^RT + cR^ Xk R^ Xk + ...) (2) 

where a, b, and c are of order Mp 2 . In this work we will only consider a spatially flat 
metric, 

d s 2 = -dt 2 + a{t) 2 5 ij dx i dx i . (3) 

The field equations are 

G>„ = T„ v + Ag^ (4) 

where 8ttG has been absorbed into and A. is a power series in derivatives, 
where the two derivative terms make up the Einstein tensor. The matter described by 
T M „ has the equation of state p(t) = wp{t). We introduce a dimensionless Planckian 
time 

t= J=— *M P1 t, A(t) = 2ln(a(t)), (5) 
V 6a + 2b + 2c 

and in the following we drop the tilde. For the above metric, 6a + 2b + 2c is the only 
combination of the four-derivative parameters that appear in the field equations, and 
we must assume that it is positive. We absorb a further factor of 6a + 2b + 2c into 
the terms in (fj|), thus making T^ v and A dimensionless. 

In terms of this dimensionless time we will be considering oscillatory behavior for 
A(t) with roughly unit frequencies and very small amplitude e. Thus d n A(t)/dt n m e 
and we consider a small amplitude expansion in powers of e. The linear terms are 
contained in Ga. 

G 00 [A(t)} = 0(e 2 ) (6) 
Gu[A(t)] = -£ c 2n d 2n A(t) + 0(e 2 ) (7) 

n>l 

The definitions above yield — C4 = 1. At 6 derivatives and above these linear terms 
originate from terms in the action containing two curvature tensors, or contractions 
thereof, and the appropriate number of covariant derivatives. 

To describe our results at higher orders in e we make the following definitions, 

G 00 [A(t)}^Y,^ n G? n) (8) 

n 

Gu[A(t)] EE -J2 C 2nQ ( p n) (9) 
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where G^ and G^ contain 2n derivatives. 

We display in the Appendix the results up to 12 derivatives, for Gp 2n ^ up to 0(e 2 ) 
and Gf ,n " > up to 0(e 3 ). As the number of derivatives grows the number of possible 
terms in the action that contribute grows very fast. Despite this we see that in 
addition to c<m only one new parameter \2n appears at each order in the total number 
of derivatives, to the order in e displayed. (Actually at 2 and 4 derivatives there are no 
parameters given the definitions above and the terms at all orders in e are displayed.) 

These results were obtained without obtaining the general field equations for gen- 
eral metric. Since the t dependence of the metric (|^) appears as a general function 
of t in g u , we can insert the metric (|3p into the action and vary with respect to A(t). 
This will give the ii field equation, and thus the Gp 2n \ for this particular metric. The 
g(? n ) can then be obtained through the use of the conservation equation (generalized 
Bianchi identity) D^G^ V = 0, which we satisfy to 0(e 3 ). In fact this completely 
determines Gp 2n ^ at order 0(e m+1 ) given Gj? n ' at order 0(e m ). All the terms in the 
field equations, to any order in e and for any given action, can be obtained in this 
way, but only the terms at the orders in e relevant for our study are given in the 
Appendix. 



2 Pure oscillations 



We start with the oscillating vacuum solutions of the linearized equations, corre- 
sponding to the Lagrangian 



C = -~A(t)(J2c 2n d?)A(t). 

n>l 

For A(t) of the form cos(u;i) the allowed frequencies must satisfy 



E (-!) nC 2n^ 



2n 



0. 



(10) 



(11) 



n>l 



When a negative or complex u 2 solves (|ITD the theory is exponentially unstable, a 
case we will not consider further. For all roots, besides the vanishing one, to be real 
and positive it is necessary (but not sufficient) that all C2 n > 0. 

Given the set of allowed frequencies cjj we can write the Lagrangian as 



C = — 2 A(t) 



i>0 



A(t), 



(12) 



after absorbing a constant normalization into A(t). We will assume that ujj > Ui for 
j > i, and uq = 0. This can be rewritten as a sum of conventional free field terms by 
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defining]] 

(hit) 



2\ 



U(d 2 + ^ 



A(t) 

(13) 



11 («? - < 
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Note that A(t) ~ cos(u>it) corresponds to all <f>j(t) vanishing except for </>i(t). The 
Lagrangian then becomes 

c = -W{-) l Ht){d 2 t +^)Ht)- (14) 

z i>0 

The (— ) l factor means that every other term comes in with the wrong sign, and in 
particular this is the case for the mode with the lowest nonvanishing frequency U\. 
The wrong sign in the Lagrangian means that the classical energy of these modes is 
negative. The energy density of the ith oscillating mode with amplitude £j is thus 
proportional to (— ) l sf. We will refer to this as a gravitational energy. These minus 
signs appearing in the linearized theory cause problems in the attempt to quantize 
higher derivative gravity. 

These vacuum oscillating solutions of the linearized theory are of course not so- 
lutions of the full nonlinear theory. In particular an oscillating A(t) will generate 
non-oscillating terms in the field equations at 0(e 2 ), and a source must con- 
tribute at this order for a solution to exist. We will include both a cosmo logical 
constant A and matter with equation of state p(t) = wp(t). With A(t) in the form 

A(t) = ecos(cj t) + a 2 e 2 cos(2 out) + ... (15) 

we expand the field equations to 0(e 2 ) using the results for and in the 
Appendix. A solution would consist in relating e and a 2 to the gravitational sources. 
Solutions at higher order in e could be found by continuing the expansion of A(t) in 
the obvious way.0 

Despite the complexity of the expressions in the Appendix, we find very simple 
results for the oscillating A(t), and a clear pattern emerges in the terms with 2 through 
to 12 derivatives. We assume that this pattern will continue at yet higher derivatives. 

e 2 -(B + Acos(2ut)) = p{t) + A (16) 
8 

9 3 

-eAcos(ujt) + e {-A + Ccos(2cut)) = p(t) - A (17) 



1 This generalizes the definitions in j| . 

2 Solutions at 0(e 3 ) were found in the 2 + 4 derivative case in 
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A = ^(-) n c 2 „u; 2 



n>l 

^=E(-) n_1 ( 2 ^- 1 ) c 2^ 2n (19) 

n>l 

C = C(a 2 ,c 2n ,X2n,u) (20) 

We see that pit) at 0(£ 2 ) is completely determined by the same constants Cj that 
appear at linear order; the additional \2n parameters have all dropped out of this 
equation. The terms linear in e in the p(t) equation must vanish separately, thus 
determining the allowed frequencies as before: uj G {uji} such that ^4 = 0. Thus at 
0{e 2 ) we find that p(t) is a constant. 

In the p(t) equation only the C cos(2u;t) term remains, which depends linearly on 
a,2 and the \2n parameters. The solution, to be compatible with p = wp, is thus 
completed by choosing a 2 so that C = 0. In this way the X2n parameters affect A(t) 
only at order 0(e 2 ), through the value of a 2 . The result is that the total pressure, 
p(t) — A, vanishes and thus 

pit) = A/w. (21) 
Then to obtain a solution for the mode with frequency cjj, its amplitude E{ is deter- 



mined from (|T6|) as 



e 



2 3 E(-) n_1 ( 2 ^ " i)^, 2 " = P(l + w). (22) 



8 n>l 



We need the value of B = J2{~) n l {^ n — l)c2n^ 2n for the various roots of A = 
X)(— ) n c 2n co ,2n . It turns out that B also has a set of positive roots that are interlaced 
with those of A. The result is that the sign of B alternates at each successive root of 
A, with the implication (assuming w > — 1) that 

sign(p) = (-)* (23) 

for the mode with frequency Ui. Thus the same modes that have negative gravitational 
energy also require a negative matter energy for their existence. In other words for 
these modes T 00 + Ag 00 < whereas the modes with positive gravitational energy have 
T o + Agoo > 0. In both cases Ta + Aga = 0. Thus in order not to violate positive 
energy conditions for matter, we may conclude that the set of oscillating solutions to 
the full theory does not include the modes with negative gravitational energy. 

Note that the vanishing of p(t) — A with A ^ requires w^0. On the other hand 
if A = then the same analysis would yield solutions for w — 0. Thus pressureless 
matter can act as a source of the oscillations, with the amplitude of the oscillations 
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determined by the energy density. In this case p can be freely chosen, unlike the case 
with A^O. 

Energy conservation for matter implies that p(t) must oscillate in an oscillating 
background, since p(t)a(t) 3 ^ 1+w ^ (where a(t) = e A ^l 2 is the scale factor) must be a 
constant. Thus an oscillating component of p(t) will occur at 0(e 3 ). In summary an 
oscillation amplitude e of A(t) corresponds to a source of 0(e 2 ) having an oscillating 
component of 0(e 3 ), a situation very different than the pure linearized theory. The 
implications of any of this for the development of a quantum theory will not be 
pursued here. 

3 Oscillations and expansion 

We turn now to consider the dynamics of these nonlinear oscillations in the context 
of the expanding universe. We will find that modes that were prohibited in the pure 
oscillating case become involved in a weak instability in the expanding case. This 
dynamics depends on the equation of state. 

We will have to explore this numerically in theories truncated to various finite 
numbers of derivatives. We use the equation of state to eliminate pit) and pit) 
and obtain an equation for A(t) alone. A single equation can be dealt with quite 
successfully numerically, even though it is nonlinear with high derivatives. Once A(t) 
is obtained then the other equation determines pit). A good check on the accuracy 
of the numerical integration is provided by the quantity p(t)a(t) 3 ^ 1+w \ which should 
remain constant. We have confirmed this for all the cases studied. 

We begin by truncating the equations at 4 derivatives and setting A = 0. We 
initially consider radiation domination w = 1/3. Here it is easy enough to keep all 
the terms, to all orders in e. A range of initial conditions gives rise to the positive 
energy expanding solutions of interest to cosmology. In Fig. ([]]) we plot a(t)/a(t)FRw 
and p(£)/p(£)frw f°r different ranges of t. For different initial conditions the size 
of the oscillating effect changes and the departure from the FRW solution can in 
particular be made much smaller, but the qualitative form of the plots remains the 
same. We will discuss the initial conditions further below; but for now we see that 
a(t) increases faster, and pit) falls faster, than in the FRW case. 

We can consider the departure from the FRW equation 3H 2 = p as an indica- 
tion of the importance of the higher derivative terms. This is seen in the plot of 
if 2 (t)/if (t) FRW and p(£)/p(£)frw m Fig- (@)- We see the relative smoothness of pit) 
and the growing oscillations in H 2 (t)/H(t)p RW , which implies the growing relative 



7 




8 



°300 305 310 315 320 1005 1010 1015 1020 

Figure 2: These two figures show the relative growth of the oscillations. 

importance of the higher derivative terms. It is in this sense that the oscillations are 
growing and signaling a dynamical instability. We find that the dynamics is such 
that departure of p(£)/p(£)frw from unity grows like a power of ln(t), and thus the 
instability is weak and physically interesting. 

Since the oscillations will not be directly observed it makes sense to consider 
quantities that are suitably time averaged, so as to smear out the oscillations. For 
example the apparent critical energy density of the universe will be given by p{t) crit = 
3(H(t)) 2 . It is then of interest to consider the quantity (p{t)) /p{t) ciit , and we display 
this along with (H (t)) / H (t) FRW in Fig. (|]). We therefore effectively have p/p CT i t < 1 
in a flat universe, and thus the oscillations are playing the role of dark energy. But 
an accelerated expansion is not a result. 

We now consider matter domination w — 0. Here we find that the oscillations 
persist, but their relative importance does not grow. For example the size of the oscil- 
lations in the quantity H{t) 2 / H{t) Fim does not grow. This is reflected by the leveling 
off of p(t) I /?(£)frw & t some value determined by the size of the oscillations. Both the 
initial and final value of p(t)/p(t)FRW is adjustable through the initial conditions. At 
the same time (H(t)) /H(t) FRW tends to unity. The physically interesting quantity 
(pit)) I p{t) crit displays the least time dependence, and its essentially constant value, 
less than unity, is determined solely by the size of the oscillations. Thus oscillations 
that have grown during a radiation dominated era become "frozen in", in relative 
importance, during the matter dominated era. 

Pressureless matter w = turns out to be the dividing line between the growth 
of the oscillations, in relative terms, for w > and the fading away of the oscillations 
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Figure 3: The bottom curve shows the smeared matter energy density (p{t)) compared 
with the apparent critical value for a flat universe, where the latter is determined by 
(H(t)). 



for w < 0. In particular a positive cosmological constant eventually dominates the 
cosmological evolution and drives the effective w negative. When this happens the 
oscillations fade away. 

How natural it is to have the oscillations depends on how close to the Planck time 
the initial conditions are specified. The initial conditions may be characterized, for 
example, by how much the acceleration parameter deviates from the FRW value at 
that initial time. A given fractional deviation specified at one initial time produces 
much less oscillation than the same fractional deviation specified at an earlier initial 
time. Indeed at an initial time close to the Planck time tpi, significant fine tuning is 
required to produce a solution close to the FRW one. 

More appealing is to consider an initial time of order 30tpi. This was used to 
produce the w = 1/3 results in Figs. where a large fractional deviation, a 

factor of 10, in the acceleration parameter was used. For a small fractional deviation, 
say at the 10% level, the resulting oscillations are tiny, and they may not grow to 
be significant until cosmologically interesting times have elapsed. This is based on 
numerical analysis that indicates that the quantity (p(t)) /p(t) ciit in the radiation 
dominated era can be parameterized as 1/(1 + cln(t) a ) for small c depending on 
the initial conditions and a ~ 5. This is for times accessible to numerical analysis, 
whereas the era of radiation domination ends at t ~ 10 28 tpi, and so whether this 
approximation remains true throughout this era is uncertain. Clearly the main issue 
is to use the time of nucleosynthesis t ~ 10 18 tpi to put a constraint on the deviations 
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Figure 4: Contributions to p(t) / p(i)FRW from terms with given numbers of derivatives. 
The roughness of the curves is an artifact of the smearing. 

from standard FRW evolution. 

We have studied the higher derivative equations in a similar way, using the results 
of the Appendix. This maintains energy conservation to 0(e 3 ), and we have verified 
the insensitivity to the X2n parameters. We will not display these numerical results, 
because they resemble the results from the four derivative equations. The new quali- 
tative feature is the existence of more than one frequency mode of oscillation, where 
the relative size of the coefficients C2 n determine how easy it is to excite the various 
modes simultaneously. The resulting oscillations can deviate quite dramatically from 
being pure sinusoidal or even periodic. But for the observable smeared out quantities, 
the results closely resemble those above. 

We now look a little more closely at the origin of the result that p(t) < p(£)frw- 
In the pure oscillation case the modes with negative gravitational energy required 
negative p(t). Here it is clear that these modes are being excited and that they are 
responsible for driving p(t) < p(£)frw- in fact i n the 4 derivative case there is only 
one oscillating mode and it has negative gravitational energy. At higher derivatives 
(H|) suggests that contributions to p(t) alternate in sign depending on the order of 
derivatives. This is displayed in Fig. @ in the case of the equations with up to 6 
derivatives and cq = 1/5. The relative sizes of the contributions are consistent with 
(f22|) given that the lowest mode with u = u\ dominates. The conclusion is that 
we are linking a dynamical instability in the expanding system to the growth of a 
negative component of the total gravitational energy. 

As a side issue we mention that when the cosmological constant is negative there 
appear to be solutions with oscillations on both large and small scales 0, where 
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the frequency of the large scale oscillation is proportional to y/X. Our results now 
allow us to determine the sign of p needed to support these solutions at any order 
in the number of derivatives, and we find that it is typically negative. This is not 
surprising since the amplitude of the large time scale oscillations can be adjusted to 
zero with appropriate initial conditions; this returns us to the pure oscillation case 
we have treated above where we found negative p(t) for some modes, in particular for 
the lowest frequency. This leaves open the possibility that with appropriate initial 
conditions it may be possible to excite the modes with positive energy in such a way 
as to have oscillations on large as well as small time scales. The practical interest in 
this is unclear. 

In this paper we have been concerned with a dynamical instability of standard 
FRW expanding universes, where the matter energy density remains positive, but is 
reduced in association with the excitation of gravitation modes with negative energy. 
This is found to be a generic feature of higher derivative gravity with any number of 
derivatives. The phenomenological implications of this will have to be explored else- 
where, but the instability appears to be of physical interest due to its slow logarithmic 
growth. This instability is incompatible with a positive cosmo logical constant. This 
observation may be interest in a theory in which the effective 4 dimensional cosmo- 
logical constant is a derived quantity, perhaps emerging from a higher dimensional 
theory || . The instability may provide a dynamical reason for the 4D flat spacetime 
to be preferred. At the same time the oscillations and their gravitational coupling 
to matter may provide a mechanism by which positive vacuum energy is converted 
through particle production to matter energy. 

Appendix 

We list here results for the quantities defined in to the required orders in the 

small amplitude expansion, and up to twelve derivatives. 

s/ 2) = mm 2 

<s P w = (& A (t)) + !(|U(t)) 2 + 3(fA(t))(|U(t)) + !(|U(t))(f A(t)) 2 
^ (4) = -U&m) 2 + Uim)(&A(t)) + 1 {*A{m§- t A(t)y 

+ f(f^))(|U(t)) 
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£/ 6) = -AwMt)) 2 - UwMt))(&A(t)) + l(f t A(t))(§A(t)) 

+ Hw A (t))(f t Mt)) 2 - (-9 - §x 6 )(f ^(t))(|U(t))(|U(t)) 

-(3 + |x 6 )(^(t)) 3 

P {8) = {&A(t)) + (f + 6 X8 )(|U(t)) 2 + (36 + 9x*){&A(t))(&Mt)) 
+ (19 + *Xs)(&A(t))(&Mt)) + Q(w A (t))(iA(t)) 

Se (8) = -KwMt)) 2 + l(^A(t))(^A(t)) - §(|U(i))(|U(0) 
Uw A (tMA(t)) + f (f A(t)) 2 (|U(t)) - (f + ixs)(|U(t)) 2 (fU(t)) 

-(-f-lx8)(|A(0)(ftA(t))(ft^)) 
"(-f -9xs)(^(t))(^(t))(^(t)) 

o (10) = (gU(O) + (^r + Xio)(|U(t)) 2 + + §Xio)(|U(0)(jUW) 
+ (59 + lxio)(^A(t))(grA(t)) + (26 + ±Xio)(|U(0)( JU(*)) 

+ ¥(f^(0)(JU(*)) 

^ (10) = i(|U(*)) 2 - f + §(£^))(g^(0) 

- §(|U(t))(|U(t)) + §(f ^))(IU(t)) 
-(-f - §xio)(|U(0)(f ^)(IU(0) 

- (-63 - ^xio)(^(0)(^W)(I^(0) 

- (-Y ~ ToXio)(^A(t))(^A(t)) 2 
- (-? - ^Xio)(SA(t))(^))(^(0) + 9(^A(t))(|A(t)) 2 
(33 + §Xio)(lU(t)) 2 (|U(t)) - (f + ± X n)(&Am&Mt))(&Mt)) 

-(f + ^Xio)(|U(*)) 2 (|U(*)) 



(12) = (&A(t)) + (318 + x 12 )(|U(t)) 2 + (549 + § Xl2 )( * A(f))(|^4(f)) 

+ (T + gXi2)(feA(t))(^(t)) + (¥ + ?xi 2 )(|U(*))(|U(0) 
+ (48 + |xi 2 )(|U(i))(j£U(*)) + 9(|A(0)(j£U(0) 

^ (12) = -f(|U(*)) 2 + §(|U(*))(IU(*)) - §(|U(*))(|U(i)) 
+ §(|U(*))(|U(i)) - §(|^(0)(|U(*)) + Uf t A(t))(£rA(t)) 
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- C-f + |x 12 )(|U(t)) 2 (|U(t)) - (¥ + ixi 2 )(|U(t)) 2 (|U(t)) 

-(-¥-|xi 2 )(fe^(t)) 2 (^(t)) 

- ~ ToXi2)(iA(t))(-^A(t))(^A(t)) 
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